Dipole-quadrupole interactions and the nature of phase III of compressed hydrogen 
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A new class of strongly infrared (IR) active structures is identified for phase III of compressed 
molecular H2 by constant-pressure ab initio molecular dynamics and density-functional perturbation 
calculations. These are planar quadrupolar structures obtained as a distortion of low-pressure 
quadrupolar phases, after they become unstable at about 150 GPa due to a zone-boundary soft 
phonon. The nature of the II-III transition and the origin of the IR activity are rationalized by 
means of simple electrostatics, as the onset of a stabilizing dipole-quadrupole interaction. 



The nature of the high-pressure phases of solid hy- 
drogen is of long time interest J5|]. One of the most in- 
triguing and outstanding issues is the nature of the 150 
GPa phase transition || and the origin of the strong 
infrared vibronic activity observed in the high pressure 
phase [[| . The present view of the low temperature phase 
diagram of hydrogen is the following. In the gas and di- 
luted phases the homopolar H2 molecules interact via 
electric quadrupole-quadrupole interactions (EQQ). At 
low temperatures, quantum rotational effects overcome 
librational barriers, so that the molecules rotate freely. 
Upon compression, EQQ interactions freeze the molec- 
ular rotations into an ordered broken-symmetry phase 
(BSP, or phase II). Above ~ 150 GPa, a third phase 
(H-A, or phase III) is observed. Optical measurements 
across the II-III phase boundary indicate a substantial 
increase in the vibron infrared (IR) activity, the effective 
charge being one order of magnitude larger in phase III 
than in phase II Q . The lack of a Drude edge in extended 
IR studies down to 1000 cm -1 indicates that phase III 
is still insulating ||. More recent optical experiments 
focused on the low-frequency region of the vibrational 
spectrum revealed the sudden appearance in phase 
III of at least twelve modes below 1000 cm -1 , whose fre- 
quencies strongly increase with pressure, extrapolating to 
zero at P < 150 GPa, and suggesting that quantum li- 
brational effects are frozen out and are not of paramount 
importance in phase III Q. 

A satisfactory theory of the origin of the IR activity 
and the nature of the II-III phase transition is still un- 
available. First-principles calculations, although numer- 
ous, have failed so far to provide a unique description of 
phase III 8 13 , even after inclusion of proton quantum 
Some of them failed in properly sam- 



effects [14 
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pling the electronic Brillouin zone [|T^,^5[|6), others in 
considering unit cells too small to be consistent with the 
multiplicity of low frequency modes |||)|,[ll]-[l4| . In par- 
ticular, some structures studied were metallic fliq] , and 
the Cmc2i structure, which is currently being considered 
as a promising candidate for phase III |) 13 TtJ, turns out 
to be unstable against Pca2\, at least within LDA and 



GGA. 

In this letter we report the results of a combined struc- 
tural and vibrational study using constant-pressure ab 
initio molecular dynamics simulations (CPAIMD) |l§| ] 
and density functional perturbation theory (DFPT) cal- 
culations |L9| . These show that, starting with the likeliest 
structure proposed for phase II, a zone-boundary phonon 
instability sets in at about 150 GPa. The resulting stable 
structures have larger unit cells and exhibit enhanced IR 
activity, precisely the features typical of phase III. We 
then show that the transition from low to high IR activ- 
ity can be understood in terms of simple electrostatics, 
as a structural transition driven by dipole-quadrupole in- 
teractions. 

The details of the ab initio calculations are the follow- 
ing. Electronic correlations were treated within a gradi- 
ent corrected density functional p0| , while the proton- 
electron interaction was described through a local pseu- 
dopotential that requires an energy cutoff of 60 Ry 
for the plane wave expansion. Convergence on BZ sums 
required 512 fc-points in the full BZ |2^] of 4 and 6- 
molecule supercells, and 256 fc-points in an 8-molecule su- 
percell. In the MD simulations, full BZ sampling was im- 
plemented using a self-consistent k-p technique |l^,|23|,Q . 

Previous CPAIMD simulations indicate that the 
quadrupolar Pca2\, hep structure (see Fig. 1) is likely 
to correspond to the low-pressure phase II Q , confirm- 
ing that the EQQ are dominant at low densities. In 
the present work we have computed the zone-center and 
zone-boundary ab initio phonon frequencies at 140 GPa 
using DFPT, and found that this structure is, in fact, 
stable. However, increasing pressure to 180 GPa, the 
phonon dispersion of the Pca2\ structure displays an in- 
stability at the (100) Brillouin zone boundary |p5|j . Us- 
ing the corresponding eigenvector we distorted the origi- 
nal structure and minimized the enthalpy by optimizing 
lattice parameters and molecular orientations. In the 8- 
molecule structure derived from the Pca2\, and shown 
in Fig. 1, we found a tendency of the molecules to lie in 
parallel, (001) hep planes. This tendency becomes more 
pronounced the higher the pressure. 
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Once molecules lie in, or close to planes, it is natural 
to surmise that another possible candidate would be a 
structure that optimizes the quadrupolar interactions in 
a triangular, two-dimensional system. Within the planes, 
neighboring molecules are oriented at 60 degrees from 
each other. This minimizes frustration of EQQ interac- 
tions, an angle of 90 degrees being impossible in a tri- 
angular lattice. Thus, we also optimized the structure 
of the triangular structure by means of CPAIMD, in a 
6-molecule cell (see Fig. 1), and found it to be stable, 
with an enthalpy lower than that of the distorted Pca2\. 
In this way we obtained a new class of structures which 
are characterized by unit cells containing more than 4 
molecules. 

In order to make a direct comparison with experimen- 
tal data, we report in Fig. 2 the vibrational spectrum 
at zone center of the undistorted Pca2\ structure at 140 
GPa and, at higher pressures, where this structure be- 
comes unstable, the spectrum of the centrosymmetric 
triangular structure. We first notice that the discon- 
tinuous frequency drop upon transformation from the 
Pca2\ to the triangular structure, compares favorably 
with an experimental jump of about 100 cm -1 . This 
jump is connected with an intramolecular bond softenin- 
ing signalled by a sudden, w 2% increase in the molecular 
bond lengths. Actually, in the triangular structure, the 
bond length is 1.41 Afor two of the six molecules in the 
unit cell, and 1.39 Afor the other four. In the Pca2\ 
structure the bond length is 1.37 A. Concomitantly, we 
find a k, 1% discontinuous volume decrease at the tran- 
sition. In fact, the main energy decrease at the tran- 
sition is in the interproton repulsion, counteracted by 
a nearly equal increase of electronic energy. The lower 
part of the experimental IR and Raman spectra of phase 
III shows a number of modes whose frequencies extrapo- 
late to zero for decreasing pressure, between 100 and 150 
GPa, but are totally absent in phase II ||. Since op- 
tics detects only zone-center modes, these data strongly 
suggest a zone-boundary libron instability of phase II, 
with formation of a larger supercell and soft mode re- 
folding at k = in phase III. This is precisely what our 
first-principles calculations predict. The mode- hardening 
with pressure which we calculate for the lowest frequency 
mode is 4l.5^/(P - 142) cm" 1 (see Fig. 2). 

The DFPT approach also allows us to compute IR 
activities [p6f , which we report in Table I, as a func- 
tion of pressure. At low pressures three of the vibrons 
are nominally IR active, but the oscillator strengths are 
very small, with effective charges Z* ps 0.02 e. After 
the transition to the high-pressure phase, the IR activ- 
ity increases suddenly by one order of magnitude, i.e. 
Z* w 0.2 e. This large jump appears to account quan- 
titatively for the experimental onset of strong IR activ- 
ity in phase III Q). The IR activity is further shown, 
using effective charges calculated at 220 GPa, to grow 
steadily by about 1.5 x 10~ 4 e/GPa with increasing pres- 



sure, to be compared with an experimental increase of 
2.07 x 10" 4 e/GPa. Finally, there is an IR phonon mode 
at about 1600 cm" 1 , which may explain that seen in one 
set of experiments 0j. 

An intriguing fact is that we find three IR active vi- 
brons against only one observed experimentally. We note 
here that this threefold multiplicity is of strictly inter- 
planar nature. Interplanar interactions cause, e.g. the 
splitting between modes A (IR active) and B (Raman 
active), as indicated in Fig. 2. If angular correlations 
between planes were to be reduced, or eliminated, for ex- 
ample by librational fluctuations, then the AB doublet 
would collapse, and the IR activity of mode A would be 
supressed. The alternative possibility of a smaller unit 
cell (3 molecules) is not supported by our calculations, 
and seems also incompatible with the large number of 
experimental librational modes. 

We computed, in addition, the static dielectric con- 
stant, which is experimentally accessible via refractive 
index measurements 27 . We predict a discontinuous 



increase of about 10% at the 150 GPa phase transition, 
which is at variance with the decrease calculated for other 
simpler structures on the basis of a dielectric catastro- 
phe mechanism, driven by diverging dipole-dipole inter- 
actions M. 

We now address the physics underlying the II to 
III transition, which appears to be driven by dipolc- 
quadrupole interactions. In fact, a simple electrostatic 
model consisting of interacting quadrupoles and corre- 
sponding induced dipoles in a hep lattice qualitatively 
accounts for the picture emerging from the above first- 
principles results. This model is an extension of the orig- 
inal model of interacting pure quadrupoles |2§| ] with the 
addition that the H2 molecules are polarizable, and de- 
velop a dipole moment when placed in an external field 
such as that resulting from all the quadrupoles of the 
neighboring molecules. Once these induced dipole mo- 
ments set in, dipole-quadrupole (EDQ) and dipole-dipole 
(EDD) interactions add to the basic EQQ interactions. 
These depend differently on the intermolecular distance 
(1/r 5 for EQQ, 1/r 4 for EDQ, and 1/r 3 for EDD), thus 
originating a competition between different orientational 
patterns p9). 

In Fig. 3 we show the energy diagram and the be- 
haviour of the dipole moment as a function of the lattice 
constant. When the molecules are far apart the EQQ 
interaction is dominant and the Miyagi-Nakamura result 
is recovered p8[ : in the hep structure for the molecular 
centers, a Pca2i structure is favored, the molecules ar- 
ranging themselves as close to perpendicular as they can. 
As the distance between molecular centers decreases, all 
individual quadrupolar electric fields on a given molecule 
increase, inducing a molecular dipole moment. The EDQ 
interaction, is however still weak in Pca2\ because in 
this structure the quadrupolar fields due to neighboring 
molecules almost cancel, and the main contribution to 
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the energy remains EQQ. With decreasing spacing, the 
increasing dipole moment would eventually drive a di- 
electric instability 0, as the electric fields created by 
the induced dipoles interact cooperatively. As any two 
molecules approach under pressure, however, the increas- 
ing EDQ interaction causes their relative orientation an- 
gle to turn, driving a k ^ libron instability and even- 
tually stabilizing a new planar triangular structure. The 
dielectric instability is therefore preempted by a struc- 
tural transition. Remarkably, the structure so obtained 
is very close to the planar quadrupolar structure with 
in-plane quadrupolar orientations just surmised from ab 
initio calculations fl3l| . 

The quadrupolar electric field at each molecular centre 
is now between 5 and 15 times larger than in the Pca2\ 
structure due to incomplete cancellation of the contri- 
butions of the neighboring molecules. This field induces 
proportionally larger dipole moments, so that the EDQ 
interaction turns out to be more than two orders of mag- 
nitude larger, and comparable in strength to the EQQ 
interaction. At the transition distance, the static molec- 
ular dipole moment jumps by a factor of about 8, closely 
paralleling that of the dynamical dipole calculated by 
DFPT, which is related to the IR activity. The jump 
of IR activity thus signals a very central aspect of the 
physics driving the II-III transition. 

In conclusion, we have presented a theory that sheds 
light on the nature of phase III of solid molecular hy- 
drogen, and explains the origin of its strong IR activ- 
ity. The quadrupolar electric field of the H2 molecules 
induces a dipole moment in the polarizable neighboring 
molecules. The ensuing dipole-quadrupole interaction, 
gives rise to a nonuniform librational instability driving a 
phase transition to a planar structure, still exhibiting in- 
plane quadrupolar ordering, but also optimizing the new 
term. The corresponding candidate structures obtained 
from first-principles calculations possess larger unit cells, 
and reproduce many measured properties, including the 
vibronic effective charges, the vibron discontinuities, the 
librational phonon spectrum, and the insulating charac- 
ter of this phase. The exact symmetry group of phase III 
and the precise reason for a single IR mode remain open 
questions, calling for further theoretical and experimen- 
tal work. 
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FIG. 1. Schematic view of the low-pressure quadrupolar 
phase Pca2i (a), and two candidates for phase III: (b) the 
distorted Pca2i, and (c) the triangular structure. Large and 
small circles indicate that the molecule lies out of the plane, 
pointing in the direction of the larger circle. Empty and 
shaded molecules correspond to different planes. 

FIG. 2. Zone-center phonon spectrum as a function of pres- 
sure for the Pca2i structure at 140 GPa (shifted to the left 
for the sake of clarity), and the triangular structure, above 
140 GPa. The pressure dependence of the IR active vibrons 
is indicated with dashed lines, and that of the soft librational 
mode with a solid line, crossing zero at P ~ 142 GPa. 

FIG. 3. Model of quadrupoles and induced dipoles: Energy 
(left panel) and dipole moment (right panel) per molecule 
as a function of the lattice constant, in an hep lattice. The 
dipole moment exhibits a discontinuity of nearly one order of 
magnitude at the transition distance. 

TABLE I. Vibron effective charges as a function of pres- 
sure. Values are reported for the Pca2\ structure at 140 GPa, 
and for the triangular structure at 180 and 220 GPa. Exper- 
imental values were taken from Ref. [4]. 
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